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ABSTRACT 

We present near-infrared J, H , and K s images and iT-band spectroscopy of two newly discovered 
stellar clusters at different stages of evolution. Our spectra suggest the presence of massive YSOs in 
the heavily embedded cluster in the star-forming region near radio source G353.4-0.4 and an 05-06V 
star in the cluster near radio source G305+00.2. We determine a if-band luminosity function (KLF) 
for both clusters and an initial mass function (IMF) for the cluster near G305+00.2. The derived 
IMF slope is V — —1.5 if the KLF is used to derive the IMF and is T — —0.98 if the color-magnitude 
diagram and spectra are used. The more reliable CMD-based slope is flatter than the Salpeter value 
usually found for stellar clusters. We find that using the KLF alone to derive an IMF is likely to 
produce an overly steep slope in stellar clusters subject to variable extinction. 

Subject headings: open clusters and associations: general — stars: formation — stars: luminosity 
function, mass function 



1. INTRODUCTION 

Despite their intrinsic rarity and short lifetimes, mas- 
sive stars are extremely important in the evolution of 
galaxies. They play an important role in determining 
the course of the formation of less massive stars, though 
the nature of this role is still uncertain, and their stel- 
lar winds and eventual supernovae shape the interstellar 
medium. They produce most of the heavy elements in the 
universe, as well as much of the UV radiation in galaxies. 
Their rarity, combined with the effects of large Galactic 
extinctions, often results in the availability of more com- 
prehensive studies of massive stars in external galaxies, 
where the entire stellar population can be observed at 
once, than within our own where massive stars must be 
studied individually and the census of massive stars is 
still very incomplete. High optical extinction within the 
galactic plane (Ay > 20) has limited optical studies of 
mass ive stars to re latively nearby regions (R so iar ^ 3.0 
kpc. iMassevI 12003). Even within that radius, optically 
selected catalogs of O stars have been found to be in- 
complete, es pecially in star-forming regions and young 
clusters (e.g. H anson fc Conti 11995). This incomplete- 
ness necessitates the use of infrared, radio and X-ray ob- 
servations, particularly in the inner regions of the Galaxy 
and in star formation regions. The near-infrared (NIR, 
1-5 /zm) is an especially useful regime for the study of 
massive stars; the stellar atmosphere is still observed di- 
rectly, but since for example Ak — 0.11 Ay, we can ob- 
serve these stars in regions where dust, either along the 
line of sight or local to the star-forming region, makes 
them inaccessible at optical wavelengths. The discovery 
and characterization of stellar clusters observable only 
in the infrared can significantly enhance our understand- 
ing of obscured Galactic regions which harbor embedded 
massive stars or massive protostars. 

1 Steward Observatory, University of Arizona, 933 N. Cherry 
Ave., Tucson, AZ 85721 

2 SETI Institute, 515 N. Whisman Road, Mountain View, CA 
94043 

3 School of Physics, University of New South Wales, Sydney 

2052, NSW, Austr alia 

Electronic address: aleistra@as.arizona.edu, acotera@seti.org, jliebert 1 



Recent studies indicate that clusters may account for 
70-90% of star formation and that embedded clusters 
(those still partially or fully enshrouded in their natal 
molecular cloud) may exceed the numb er of more tradi- 
tional open clusters by a factor of ^20 (jElmeereen et alJ 
l2000HLada fc Ladall2003h . In the last decade, advance- 
ments in NIR observational capabilities resulted in the 
discovery and classification of some of the most mas- 
sive young stellar clusters in the G alaxy, each contain- 
ing dozens of O a nd WR stars (e.g. iNagata et al.lll995l: 
iCotera et aDll996t iFiger. Morris, fc McLeanlll996j) . Re- 
cent studies l)Figer et all I1999J) have suggested that 
within these clusters, the initial mass function (IMF) 
does not follow the canonical Salpeter form with a slope 
T = —1.35, but instead is more heavily weighted to- 
ward massive s tars: mass segregat ion has been proposed 
as a solution ( St olte et al. 1 12002) . In the last several 
years a number of studies of well-known star forma- 
ti on regions have also been carried out in the NIR, (e.g 

E kumura et ap 120001: iBlum. Damineli fc Contil 120011: 
pnti fc BlunTl2002t iFigueredo et alJ l2002j) . These 
studies have in most cases found an IMF consistent 
with the Salpeter value, and have uncovered candi- 
date massive YSOs. In addition, within the past 
ten years, massive YSOs with in molecular clouds have 
been studied in t he NIR, (e.g. iChakrabortv et alll2000t 
Ishii et all 120011) and in young stellar clusters (e.g. 
Hans on. Havworth. fc Contilll997[) . Massive YSOs, how- 
ever, remain significantly less studied and are poorly un- 
derstood in comparison with their lower-mass counter- 
parts; many more must be identified and studied before 
we can adequately address how the formation of massive 
stars differs from that of low-mass stars. 

The final release of the Two Micron All Sky Sur- 
vey (2MASS) has fostered studies which can probe the 
entire Galaxy for previously unknown stellar clusters. 
Initial attempts were made which searched for stellar 
density enhancem ents, (e.g. IDutra fc Bical 120001 120011: 
IDutra et al.l 120 031 . but the identification of previously 
unkno wn clusters has met w ith limited success. For ex- 
ample, |Dn^^^^Bica (2009) identifie d 52 candidate clus- 
'aierfeowitrMti, gi^Miph^xitnsjbgdiivait ions IjDutra et alj f2003) 
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indicated were in fact 10 confirmed clusters, 3 "prob- 
able" clusters, and 11 "dissolving cluster candidates"; 
the remainder wer e not clusters. Ou r observations of 
at least one of the iDutra et al.l (|2003T) "confirmed clus- 
ters" , however, indicates that the "cluster" is most likely 
a region of low extinct ion rather than a true cluster 
ijCotera fe Leistral 12005). We have performed an inde- 
pendent search of the 2MASS archive, using color crite- 
ria in addition to stellar density enhancements. We have 
searched in the vicinity of regions identified as likely sites 
of star formation based on radio and IRAS far-infrared 
flux ratios, and are currently conducting a search of the 
entire 2MASS Point Source Catalog. We search the Point 
Source Catalog for regions of higher stellar density than 
the background (determined locally within a 5'radius) 
which are redder in H — K than the local field. This se- 
lects for embedded clusters, with the color criteria help- 
ing to eliminate chance superpositions and regions of 
low extinction. In contrast, Dutra & Bica (2000,2001) 
use only stellar density to select clusters. Our method 
has been relatively successful to date; correctly selecting 
7 clusters out of 9 potential targets, including 4 candi- 
dates toward the inner Galaxy. We present NIR imaging 
and spectroscopy of the two confirmed clusters in the 
inner Galaxy in this pape r, and discuss the two u ncon- 
firmed targets in detail in lCotera & Leistral {2005). The 
clu ster near G305. 3+0 .2 was independently discovered 
bv IDutra et alJ (|2003b). The additional 5 outer-galaxy 
targets are described in Paper II. 

NIR imaging and spectroscopy of both young stellar 
clusters and nascent stellar clusters enables us to ex- 
pand the study of the IMF in objects where there has 
been little to no stellar evolution off the main sequence 
or cluster evaporation, and where the cluster age can 
be constrained to within ~ 2 Myr. Spectral typing of 
the most massive stars in the cluster allows their masses 
to be determined relatively precisely, and when com- 
bined with photometry it facilitates a reliable determina- 
tion of the masses of stars throughout the entir e cluster 
(Mas sey, Johnson, & DeGioia-Eastwood 1995; iMassevI 
2002) , allowing the initial mass function of the cluster to 
be determined more accurately than photometry alone 
would permit. In this paper we present the results of 
NIR observations of two clusters found toward the inner 
Galaxy, which we designate by the Galactic coordinates 
of their centers, G353.4-0.36 (17:30:28 -34:41:36 J2000) 
and G305. 3+0.2 (13:11:39.6 -62:33:13 J2000). In Paper 
II we will present the results of similar observations of 
five clusters in the outer Galaxy. 

In Sj21we present the observations and data reduction, 
in S|3 we present the spectra and classifications of the 
spectroscopically observed cluster members as well as the 
color-magnitude diagrams, and in 21 we describe the lu- 
minosity function and the initial mass function. 

2. OBSERVATIONS & DATA REDUCTION 

We observed candidate young stellar clusters with the 
facility instrument IRIS2 on the 3.9m Anglo- Australian 
Telescope (AAT) on July 12-15, 2003. IRIS2 is an 
imaging spectrometer which uses a 1024x1024 Rockwell 
HAWAII- 1 HgCdTe array with a platescale of O'.'45/pixel, 
resulting in a 7.'7x7f7 field of view. Images were obtained 
in J (1.25 jum), H (1.63 /Jim), and K s (2.14 fim) filters. 
R ~ 2300 spectra of selected stars in each cluster candi- 



date were obtained in K for all candidates. 

We selected a total of four cluster candidates in the 
southern hemisphere using the 2MASS Point Source Cat- 
alog based on color and density criteria. Two of the 
candidates observed appear to be r egions of low ex- 
tincti on and are discussed elsewhere (|Cotera fc Leistral 
2005). The two confirmed clusters are near radio H II 
regions designated G305.3+00.2 and G353.4-0.4. We 
present three-color composites of the 8' x 8' images of the 
G305.3+00.2 and G353.4-0.36 clusters in Figures d and 
respectively. 

G305.3+00.2 is an H II region which has been 
previously observed us ing radio recombination lines 
llWilson fc Mezgerll970ft . C I emission in the submillimc- 
ter (|Hua*nT^t*aLf|l9 99). and in the mid-infrared (MIR) 
by the Midcourse Space Experiment (MSX). The kine- 
matic distance of 3.5 ± 1.1 kp c obtained for this H II 
region ijWilson fc MezgeJll970(l agrees well with the dis- 
tance of 3.3 kpc fo r masers several arcminutes away 
(|Caswell et al.lll995l) . suggesting they may be part of a 
single star-formation complex. A distance of 4 kpc is 
adopted as an upper lim i t to t he radio kinematic dis- 
tance by iClark fc Portel (2004) in a study of the star 
clusters Danks 1 and 2 in this region. The situation is 
more complex for the G353.4-0.36 cluster, which is in 
a region known to be a site of massive star formation. 
There are numerous radio sources located within 1' of 
the NIR cluster, which we discuss in detail in i )3.2.2l 

All photometric observations were done in excellent 
seeing conditions: 0'.'7-0'.'9. The images were reduced 
and combined automatically at the telescope using the 
ORAC-DR pipeline. ORAC-DR is a generic data re- 
duction pipeline created at the Joint Astronomy Centre 
in Hawaii, originally for use with various UKIRT and 
JCMT instruments. Subsequent reprocessing did not no- 
ticeably improve the images, therefore the pipeline pro- 
cessed data has been used throughout. Source detection, 
PSF fitting, and photometry was carried out using IRAF- 
DAOPHOT, and is discussed in detail in $01 

All spectra were obtained with a l"x 7.'7 slit. The 
long-slit format combined with the high stellar density 
within the FOV resulted in the simultaneous observation 
of multiple stars. Total integration times ranged from 10 
minutes to 30 minutes, and were chosen to provide ad- 
equate S/N for NIR spectral classification as described 
in Hanson, Conti, & Rieke (1996). After the data was 
flat-fielded, grism curvature was removed using the FI- 
GARO tasks cdist and sdist. Wavelength calibration 
was performed using the the OH~ night sky lines and 
the FIGARO task arc. The uncertainty in the wave- 
length calibration fit was determined to be 2.18 A. The 
FIGARO task irflux was used both to flux-calibrate the 
spectra and remove the telluric absorption using the G2V 
standards HD157017 and HD115496. Both of the stan- 
dards had intrinsic Br 7 in absorption, with equivalent 
widths of 5.7 A for HD157017 and of 5.6 A for HD115496; 
in each case, the absorption line was removed by fitting 
a line to the continuum in the region of the line in the 
standard star spectrum prior to flux calibration. The in- 
dividual spectra were obtained by extracting apertures 
4-5 pixels wide from the full spectral array, then per- 

4 FIGARO is part of the Starlink software package available at 
bttp: / /star- www. rl.ac.uk/ 
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Fig. 1. — Color composite (J=blue, H=green, K=red) of the region around the G305+00.2 cluster. Image is approximately 8' on a side. 
The cluster is clearly apparent as a concentration of stars with similar colors; no nebular emission is apparent in the immediate vicinity of 
the cluster though a ridge of nebulosity is present to the northwest. 



forming background subtraction using apertures of the 
same width on either side of the source, separated by 2 
pixels (0'.'9). We also extracted off-source spectra in each 
cluster to characterize any nebular emission. 

3. ANALYSIS 

3.1. Spectroscopy 

The development of NIR spectra l atlases of nearby 
massive stars of kn own spectral type i Hanson et al .119961: 
IMorris fc flerahylllflM iRhim et a,l]ll997fl. provides a 
valuable classification scheme for stars too heavily ob- 
scured by dust to permit optical spectroscopy. In the K 
band, in addition to the Br7 (2.165 fim) line, massive O 
stars have helium (He I 2.058 ^m, He I 2.112 /zm, He II 
2.189 urn), carbon (C IV 2.078 fxm), and nitrogen (N III 



2.116 fim) lines in their spectra which allow for the de- 
termination of the spectral type to within a subtype if 
th ere is adequate (> 70) line signal to noise. Table 6 
of lHanson et alJ 1^996) indicates that in many cases the 
mere presence of these lines in emission or absorption 
(without considering equivalent width) is sufficient to de- 
termine spectral type to within two subtypes for O stars. 
The situation is more complicated for B stars, which have 
fewer features in this part of the spectrum; however, they 
are still classifiable using only if-band spectra. 

We obtained if-band spectra of five stars in the 
G305.3+0.2 cluster field and three stars in the G353.4- 
0.36 cluster. In order to reduce the level of foreground 
contamination, we imposed a color cut of H — K > 0.5 
based on the 2MASS magnitudes and selected the bright- 
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Fig. 2. — Color composite (J=blue, H=green, K=red) of the region around the G353.4-0.36 cluster. Image is approximately 8' on a 
side. The cluster is surrounded by intense nebular emission and is contained in a larger dark molecular cloud. 



est stars meeting this requirement. Despite this cutoff, 
two of the five stars observed in the G305.3+0.2 clus- 
ter proved to be foreground contaminants with sufficient 
line-of-sight extinction to push them over our threshold. 
The cluster sequence was much narrower and more well- 
separated from the foreground in the G353.4-0.36 cluster, 
and no obvious foreground contaminants were present in 
our spectroscopic sample. The G353.4-0.36 cluster was 
sufficiently red (H — K c i uster > 1.3), that the time re- 
quired to obtain a useful signal-to-noise in if-band spec- 
tra would have been prohibitively large, so only i-T-band 
were obtained. 

3.1.1. G305. 3+0.2 Cluster 



We present spectra for the three cluster members, 
which we label A1-A3, in Figure 

In Figure^Jwe present a 106" x 120" image of the clus- 
ter and label the positions of sources A1-A3. 

The measured magnitudes (see E I3.2[) and observed 
spectral lines for A1-A3 are presented in Table 2] The 
other two stars for which we obtained high S/N spectra 
have late-type spectra, as indicated by strong CO ab- 
sorption at 2.29 and 2.32 fim, suggesting they are either 
foreground objects or YSOs. The lack of nebular emis- 
sion in the cluster and the presence of weak (nearly the 
same as in the G2V spectral standard) Br 7 absorption in 
one of the spectra suggest that these are foreground ob- 
jects rather than YSOs. In addition, the K magnitudes 
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Fig. 3. — Spectra for the cluster stars in the G305. 3+0. 2 cluster. 
Top panel: Source Al, identified as 05-6V. Middle panel: Source 
A2, identified as B0-1V. Bottom panel: Source A3, identified as 
B2V-B3V. 




Fig. 4. — Region immediately surrounding the G305.3+0.2 clus- 
ter. Spectroscopically classified sources are marked as Al, A2, and 
A3 and sources showing CO absorption are labeled with "CO". 
This image is approximately 120"x 106 "and the "cluster" area is 
marked. 



of these objects (K = 9.43 and K = 10.114) make them 
too bright to be low-mass YSOs at the cluster distance, 
and the presence of main-sequence O and B stars argues 
against identifying these objects as massive YSOs. We 
thus conclude that these two stars are most likely late- 
type foreground stars, and excluded them from further 
analysis. 

Although nebular emission can be seen in the full image 
(Figure^!, it is significantly removed (>!') from the clus- 




FlG. 5. — K-hand image of the G305.3+0.2 cluster region with 
8 /im contours from the MSX mission. The K image has been 
stretched to emphasize the nebular emission. Note the close corre- 
spondence between the mid-IR emission and the nebular i^-band 
emission. 



ter. Nevertheless, in order to ensure that any measured 
Br7 (2.166 /zm) is stellar in origin and not contaminated 
by nebular emission within the cluster, we extracted a 
local background spectrum for the cluster. There were 
no features apparent in the resulting spectrum; we thus 
conclude that nebular emission within the cluster is neg- 
ligible. This conclusion is supported by an apparent bub- 
ble of MIR emission seen in the MSX Band A image (see 
Figure EJ); the MIR emission avoids the cluster itself. 

Figure 13 shows that source Al has emission lines with 
equivalent widths stronger than -2 A at 2.116 /im and 
2.166 (im (see Tabic The line at 2.166 /im is im- 
mediately identifiable as Br7. We identify the line at 
2.116/im as N III, which is consistent with the lines used 
in the the classification system presented in lHanson et alJ 
(1996); the broad nature of this line is due to the mul- 
tiplet nature of the transition responsible rather than 
broadening by stellar winds. The presence of Br7 and 
N III 2.116/im in emission, without further information 
and without equivalent widths, is sufficient to identify 
the star as being an early to middle O supergiant; the 
broad Bi'7, produced in th e stellar winds, is no t observed 
in main-sequence O stars l)Hanson et al. I ll 996'). There is 
a possible weak detection (~ 2a) of C IV in emission 
at 2.078 /im. Th is line only appears in O stars ranging 
from 05 to 06.5 ijHanson et al.ll 9961. and if real, signifi- 
cantly constrains the stellar type. Helium lines arc often 
observed both in emission and absorption in the spec- 
tra of massive stars: He I (2.058 /im), He I (2.112 /im), 
and He II (2.189 /im), are all absent from the spectrum 
of Al. Poor removal of the telluric features near the 
2.058 /tmfeature prevents us from drawing any conclu- 
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sions based on our non-detection. If real, the absence of 
the He I (2.112 /im) line restricts the spectral type to 
06 or earlier. A He II line is expected in an O star; by 
estimating the strength of possible features dominated 
by the noise (as described in detail in § 13.1.2(1 we can 
place an upper limit of 0.5 A on the equivalent width of 
any potential He II (2.188 /im) feature. This is consis- 
te nt with the width of the feature in the stars observed 
bv lHanson et al.l (J1996) , so the non-detection does not 
rule out an O star identification for this source. Taken 
together, these spectral characteristics suggest a spectral 
type of 05Ib-06Ib for Source Al. If the weak detection of 
C IV is discounted, the presence of the N III line and the 
limit on an He II line at 2.188 /im allows an 07-08 iden- 
tification as well. Even when present, however, the C IV 
line is weak, with an equivalent width weaker than -2 A; 
thus, while a positive detection of this line would allow 
for definitive classification of this source as an 05Ib-06Ib 
star, a non-detection at the given S/N does not preclude 
the same classification. 

The intrinsic NIR colors of O and B stars range from 
-0.08 to -0.01 l(WegnerllT99l : this small range allows 
an extinction to be derived even without knowing the 
precise spectral type of a massive star. For source 
Al, the extinction thus derived based on the observed 
H — K color is Ay = 12 assuming the extinction law of 
iRieke fc Lebofskvl l| 19851) . However, the large range in 
absolute Mk for O supergiants prevents us from making 
a distance determination based on Source Al. We can 
only say the distance is greater than ~ 3.3 kpc, which 
would be the d i stance for a main-sequence 05-06 star. 
iClark Ez Porter! l(2004[) adopt a distance of 4 kpc to the 
Danks 1 and 2 clusters in the same star formation com- 
plex, calling it an upper limit to the values allowed by 
the radio and Ha observations, and we will follow suit, 
acknowledging that the uncertainties in this value are 
~ 0.5 kpc. 

Source A2 shows a strong Br7 (2.166 /mi) line in ab- 
sorption with an equivalent width of 6.2 ± 1.2 A and a 
probable weak He I (2.112 /im) line in absorption with 
EW = 0.7 ± 0.2 A. This combination of features occurs 
only in B stars; a comparis on of the equivalent width of 
the lines with the B stars of lHanson et alJ l|1996f) suggests 
an spectral type in the range of B2-B4. If the He I line is 
considered only as an upper limit, the classification be- 
comes more problematic, and the star could range from 
B2-A2. The star has H - K = 0.68, which for any star 
in this range of spectral type excludes a foreground ob- 
ject. Unlike for Al, the luminosity class of these sources 
cannot be de t ermin ed from these spectral features; as 
lHanson et al.1 l(1996t) points out, the A-band spectra of 
early B supergiants are indistinguishable from those of 
early B main-sequence stars about half the time, and 
those of late-B supergiants cannot be distinguished from 
early-B dwarfs. 

If we assume that A2 is a cluster star, we can con- 
strain the absolute magnitude, and thus the spectral 
type, by requiring the distance to be the same as for 
the O star. Since the intrinsic near-infrared colors vary 
by less than 0.1 magnitude for stars i n the range of s pec- 
tral types allowed by the spectrum fWegnerl 1 1994ft . we 
can derive a extinction for this source rather than use 
that derived from the O star, thus reducing the effects 
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Fig. 6. — Spectra for sources in the G353.4-0.36 cluster. All 
three are identified as massive YSO candidates. The Br7 emission 
line seen in Bl is contaminated by nebular emission (see Fig.gJ. 



of differential extinction. This gives an extinction to 
source A2 of Ay = 11.6, or Ak = 1 -3 using the red- 
dening law of IRieke fc Lebofskvl l(1985|) . At the distance 
of 4 kpc, we obtain an absolute magnitude for Source A2 
of Mk = —4.0, roughly that expected for an 08V star. 
This identification is not consistent with the spectral fea- 
tures of A2; a smaller distance, or an identification of A2 
as an early-B supergiant, could explain the spectrum of 
A2. If the radio distance of 3.3 ± 0.3 kpc is used in- 
stead, we obtain an absolute magnitude of Mk = —3.3 
for Source A2, making it a B0V-B1V. 

Source A3 shows only Br7 in absorption with an EW of 
5.9±1.3A. We place an upper limit on an He I absorption 
line at 2.112 /im of 0.6 A. As discussed above, this width 
for Br7 only constrain the classification of the star as 
main sequence B or early A. The observed K magnitude 
is 11.96, which corresponds to an absolute Mk — —2.6 
assuming the extinction and distance of an 05Ib-06Ib 
star for source Al; this is consistent with an identification 
of A3 as a main-sequence B1V star. The radio distance 
would imply a B2V identification, also consistent with 
the spectral features of A3. Source A3 is not among the 
brightest stars in the cluster region; it happened to fall in 
the same long slit as one of the foreground contaminants 
we had targeted for observation. This suggests that the 
other cluster members brighter than A3 are also late O 
or early B stars. 

3.1.2. G353.4-0.36 Cluster 

Spectra for the three sources observed in this cluster 
are presented in Fig. 

An enlarged version of the relevant portion of Fig. [3 
is presented in Fig. with the positions of the spectro- 
scopic targets indicated with arrows and labels. 

The only non-nebular feature which we detect is CO 
absorption in Source Bl; the Br 7 emission observed in 
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Fig , 7. — Maser positi o ns from the l i teratu re ICaswell et alJ 
1200(1) : lArgon et all <2000h : IVal'tts et all 1200(3) ') overlaid on the 
G353.4-0.36 cluster K-band image. Note that they appear in re- 
gions which are dark in the near-IR, suggesting a more deeply em- 
bedded origin. Sources B1-B3 arc indicated, and all cluster sources 
detected in H and K are marked with crosses. 
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Fig. 8. — Nebular spectrum from the G353.4-0.36 cluster region. 
Emission lines present are He I 2.058 fim, H2 2.12 (im, and Br7 
2.166 [irci. 



all three spectra is contaminated by nebular emission 
to such a degree that we cannot disentangle any stellar 
component that may be present. While this line is much 
stronger in Bl than in the other two sources, the nebular 
emission is highly spatially variable in the cluster region 
and this does not demonstrate a stellar origin for the line. 
Additionally, the line width is significantly narrower than 
that of Source Al and is similar to that observed in the 
off-source nebular spectrum (Fig. |SJ). 

The CO absorption in Source Bf in combination 
with the red colors (Table |3} are similar to those as- 
sociated with solar-m ass young stellar objects (YSOs) 
ijGreene fc Ladallf 996|) . or a cool giant or supergiant. If 
Bf is a YSO, the CO absorption is from the circumstellar 
material; otherwise it is pho tospheric in nature. Us ing 
the radio kinematic distance l|Forster fc C aswell 2000) of 
3.6 kpc to the cluster, we derive an Mk for Source Bf 
of -0.8 without correcting for extinction. Correcting for 
extinction is difficult to do accurately in this region of 



highly variable extinction, especially when the intrinsic 
colors are not known since the nature of the object is un- 
certain. Nevertheless, limits can be placed on the amount 
of extinction present, and thus the absolute magnitude 
of Source Bf . The lower limit is given by the uncorrected 
value of Mk = —0.8, which assumes the color observed 
is the intrinsic color, while the bright limit can be de- 
rived assuming an intrinsic H — K = 0.3, characteristic 
of late-type stars; this gives an extinction to source Bf of 
Ay = f6.6 magnitudes and an extinction-corrected ab- 
solute Mk of —2.6. This is several magnitudes brighter 
than the expected magnitude of YSOs of approximately 
a solar mass at the distance and extinction of this cluster, 
Mk ~ 1— 3 ijOasa. Tamura. fc Sugitanill999|) . and some- 
wha t lower than the Mk for massive YSOs, Mk <~ — 1 to 
—5 l|Ishii et al.l 200f). Finally, we n ote that this Mk is 
consi stent with that for a 7M Q YSO l)Chakrabortv et all 
12000(1 . We conclude that if Source Bf is a YSO, it has a 
mass greater than a few solar masses based on its abso- 
lute magnitude in K , but observations of more massive 
YSOs are still sufficiently few that a more accurate mass 
determination based solely on the absolute magnitude is 
not possible. Given the nebular emission, seen as He I 
(2.058 fim), H2 (2.12 /xm), and Br7 (2.f66 /mi) emission 
off the stellar sources (see Figure |SJ|, G353.4-0.36 is ob- 
viously a region of current star formation; therefore, the 
identification as a massive YSO is more probable than a 
late type cool giant or supergiant located in the cluster 
itself. 

Since Source Bf was not detected in J, it cannot be 
placed on a color-color diagram to determine whether a 
NIR excess is present, which could help to discriminate 
between the YSO and cool field star possibilities. For 
Bf to be a cool giant, it would need to be a foreground 
star with the appropriate color and magnitude, which 
falls by chance in the cluster region. Rather than use 
the entire 8'x 8' field to determine the field star den- 
sity, as we did for the G305.3+0.2 cluster (JOJ, we used 
only the heavily extincted region surrounding the clus- 
ter. This is because the molecular cloud in which the 
cluster is embedded extinguishes the background stars 
to such a degree that using the entire field would signifi- 
cantly overestimate the level of field star contamination 
in the immediate region of the cluster. We estimate the 
probability of a field source as bright as or brighter than 
Source Bf and red enough to satisfy the color cut falling 
within the cluster region to be approximately 18%. This 
is a conservative estimate, since at the edges of the cloud 
reddened sources become visible and increase the field 
star density, especially of red objects, over what it would 
be at the location of the cluster. Nevertheless, we cannot 
rule out either a foreground giant or a YSO explanation 
for Source Bf . 

As with Source Bf , the non-detection of Sources B2 
and B3 in J prevents us from using a color-color diagram 
to measure NIR excess. No photospheric features are de- 
tected in the spectra of either Source B2 or B3; Source 
B2 shows a rising spectrum in K suggesting a strong 
NIR excess, while the spectrum of B3 is essentially flat 
in this region. In order to determine whether the spec- 
tra were truly featureless or merely had a signal-to-noise 
too low to see expected features, we fit a continuum to 
the spectra and examined all excursions above and below 
the fit. 90% of these deviations had an equivalent width 
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less than 1.7 A. Fo r comparison, the detec ted absorption 
lines tabulated by iGreene fc Ladal l)1996l) for low-mass 
YSOs range in equivalent width from 0.3-5.6 A for Na I 
and Ca I, with C O usually exceeding 2 A when present. 
Ilshii et all l|2001|) conducted a similar survey of massive 
YSOs; the only emission lines other than Br 7 detected 
in a significant number of sources are CO (with an equiv- 
alent width exceeding 4 A) and H2 (with an EW > 3 A 
in all cases, and > 5 A in most cases). We thus con- 
clude that Source B2 is genuinely featureless, but cannot 
classify it. The final source, Source B3, had no reliably 
detected features but the signal-to-noise was low enough 
that we cannot reliably call it featureless. 

The observed K magnitudes are consistent with a 
B star identification for sources B2 and B3; how- 
ever, the extincted but distance-corrected Mk magni- 
tudes of ~ —0.2 to —0.6 are also similar to those ob- 
served for the massive Y SO (M ~ 7M Q ) 05361+3539 
ijChakrabortv et aT1l2000D . Thus, although these sources 
are massive, we cannot distinguish based on their NIR 
spectra or magnitudes between shrouded B stars and 
less-evolved YSOs. Mid-IR observations with sufficient 
resolution to resolve the individual sources (separated 
by ~ 5") would aid in this determination; deeper J-band 
photometry, detecting more of the cluster stars, would 
also be useful. We note that although we see ionized gas 
suggesting the presence of O stars, we have not detected 
any O stars which would be the source of the ionizing 
radiation in this cluster. 

Due to the young age of the sources observed in this 
cluster and the lack of photosphcric features in their spec- 
tra, the spectra were unsuitable for determining a reliable 
distance. Thus, th e kinematic distance to t he associated 
maser and UCHII ( Fors ter fe Ca swell 2000J) was used in- 
stead, adjusted to a distance to the Galactic Center of 8 
kpc from the original 10 kpc. This gave a distance to the 
cluster of 3.6 kpc. Assuming an intrinsic H — K = 0, we 
estimate the reddening to the cluster at Ay = 22 based 
on the narrow cluster sequen ce at H — K ~ 1.3 and as - 
suming the extinction law of iRieke fc Lebofskvl |l985). 
This estimate is highly uncertain due to the young age 
of the sources; many are likely to have a near-infrared 
excess leading to an overestimate of the line-of-sight ex- 
tinction to the cluster. 

3.2. Photometry 

We obtained images in J, H, and K s of both clusters to 
a limiting magnitude of approximately J = 16, H = 18, 
K s = 18.5, with total integration times of 12 minutes in 
each band. The limiting magnitude was brighter than 
expected due to confusion, which is most noticeable in J 
due to the slightly larger PSF and the greater sensitiv- 
ity of the instrument at shorter wavelengths. Seeing was 
/ .'7-0'/8, which, since the IRIS2 platescale is 0'.'45/pixel, 
resulted in a slight undersampling of the point spread 
function (PSF), thus making PSF fitting more uncertain. 
Our individual images were taken using a random dither 
pattern with sub-pixel dithers employed to improve the 
PSF. In an effort to better understand our errors we per- 
formed both PSF fitting and aperture photometry for 
each source. There was no systematic offset between the 
two methods, but the errors were ~ 2 times larger for 
the aperture photometry due to the crowded fields. 



Photometric calibration was performed using the 
2MASS magnitudes of field stars, after correcting from 
the IRIS2 fil ter system to the 2MASS filter system as 
described in iCarpenterl l)2003j) . The calibrated magni- 
tudes for the stars in the cluster area are presented in Ta- 
ble[21 The large field of view and location in the Galactic 
Plane provided over 100 stars in each pointing which were 
bright enough to have good photometry with 2MASS, 
but faint enough to be unsaturated in our IRIS2 images 
(11.5 < K s < 14). Those stars which were relatively 
isolated in the IRIS2 images were used as the photomet- 
ric calibration set. We chose to use a relatively large 
number of calibration stars rather than selecting the few 
most isolated stars to reduce effects of potential variabil- 
ity and photometric outliers among the calibration stars. 
The scatter in the photometric calibration derived from 
comparison to 2MASS is the dominant source of photo- 
metric error, contributing two to three times the mea- 
surement errors as reported by DAOPHOT. DAOPHOT 
errors were ~ 0.03 mag while the calibration uncertain- 
ties were (A J = ±0.05, AH = ±0.06, and AK = ±0.06 
mag). Quoted errors in the 2MASS photometry were 
negligible, with most stars having an error of ±0.003 mag 
or less in all bands. Thus, the quoted error should be 
considered an overestimate when considering the relative 
photometry of stars within either cluster; the calibra- 
tion errors from comparison to the 2MASS photometry 
will shift all our measurements by the same amount. No 
trend in the photometric errors, either internally or rel- 
ative to the 2MASS data, was observed with location. 
Finally, the positions of the stars were also adjusted to 
agree with 2MASS by minimizing the offsets between the 
2MASS and IRIS2 positions allowing for pointing offset 
and rotation. 

3.2.1. G305. 3+0.2 

The color composite of the full J, H, and K s images is 
presented in Fig. ^ the cluster alone is shown in Fig. 21 
with the spectroscopic targets marked. 

The cluster is clearly visible in the full-size image with 
a concentration of nebular emission to the northwest. In 
order to help determine whether the nebular emission is 
physically associated with the cluster, we overplotted the 
contours at 8 /mi from the MSX mission 5 . (Fig. 0). The 
ridge of near-IR nebulosity corresponds to the brightest 
portion of a roughly circular structure of mid-IR emis- 
sion, with the cluster located in the interior where there 
is no mid-IR emission present. The general appearance 
is that of a wind-blown bubble, and the 8 /^m emission 
wraps entirely around the cluster at a lower level. The 
cluster is located off-center in this structure, near the 
brightest portion of the mid-IR emission, but there is 
no mid-IR emission and no near-IR nebulosity present 
in the area of the cluster itself. The cluster is dense and 
well-defined, with stellar density much higher than in the 
field. 

The K versus H — K color-magnitude diagram of the 
cluster region is shown in Figure El 

At radii of approximately 30" in the east-west direc- 
tion and 20" in the north-south direction from the cluster 
center the stellar density has fallen to that of the field, 

5 On-line data are avail able from 

http://www.ipac.caltech.edu/ipac/msx/msx.html 
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Fig. 9.— K vs H-K for the G305. 3+0.2 cluster region, with 
all sources plotted. Typical error bars are smaller than the circles. 
Filled symbols designate spectroscopic targets A1-A3. 




Fig. 10. — K vs H-K for a randomly selected control field for 
G305.3+0.2, with the same area as the cluster field. 



which we used to define the cluster region. Foreground 
stars are apparent in the color-magnitude diagram at 
H — K 0.3; in this cluster there is no clear separation in 
color between cluster and field stars, just an overdensity 
of redder stars in the cluster; as a result, we cannot im- 
pose a firm color cut to separate field stars from cluster 
stars. A color-magnitude diagram of a randomly selected 
control field with the same area as the cluster is shown 
in Fig. 1101 many fewer stars are present, especially at 
bright magnitudes and moderately red colors. 

In order to account for field star contamination within 
the cluster region, we determined the average number 
of stars per square arcminute in the image outside the 
cluster region in color-magnitude bins of AK = 0.5, 
A(H — K) = 0.5 and randomly selected the appropri- 
ate number of stars from the cluster field for removal. 
This is similar to the procedure employe d by, among oth- 
ers, B lum, Conti, & Damineli (2000) and lFigueredo et al.l 
(2002). In cases where less than one star was expected in 
the cluster field in a particular color-magnitude bin, the 
number expected was used as a probability for removing 
a star. A total of 24 "field" stars were removed, leaving 
115. The main concentration of cluster stars is at about 
H — K = 0.8, with a gradually declining number present 
out to H - K — 4. The resulting cluster CMD with the 
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Fig. 11. — Distance- and reddening-adjusted K vs. H-K for 
the G305.3+0.2 cluster region, stati stically corrected for field star 
contamination. The ZAMS from the lMevnet fc Maederl 120031) evo- 
lutionary models has been transformed to the observed quantities 
and overplotted. 



field stars statistically removed is shown in Fig. ITTI 

Given the spectroscopically confirmed presence of OB 
stars in the cluster, as well as the lack of an obvious 
color gap, we consider it more likely that these very red 
sources are either background sources or sources with a 
near-IR excess due to local dust than that they represent 
a separate cluster giant branch. The red sources are not 
concentrated toward any part of the cluster, though they 
may occur more frequently on the outskirts (as would be 
expected if they are background objects). Sources redder 
than H—K = 1.5 were excluded from analysis of the clus- 
ter KLF and IMF; they are unlikely to be main-sequence 
cluster members. If they are included and assumed to 
be on the main sequence, the resulting extinction correc- 
tion would give very large values for the masses and an 
overly flat slope to the IMF. If these sources are clus- 
ter members, they are pre-main-sequence objects, and 
their masses are difficult to determine from H and K 
photometry alone. Thus, including them in the IMF de- 
termination would give an inaccurate result whether or 
not they are cluster members, and they have been ex- 
cluded. Finally, the crowded nature of the cluster region 
means that these very red sources may suffer from poor 
photometry. 

The J — H vs. H — K color-color diagram (Fig. I12JI 
is of limited utility in identifying cluster members or de- 
termining whether some cluster members are pre-main- 
sequence objects. 

Since many sources were undetected in J, it will not 
represent all cluster members, and faint red sources 
(where we would expect to find the relatively low-mass, 
pre-main-sequence objects) would be most commonly 
missed in the color-color diagram. A cut based solely 
on H — K must still be applied to exclude background 
sources. Fig. El shows few sources in the area occu- 
pied by pre-main-sequence objects. Of those sources 
separated from reddened main-sequence stars by more 
than 3<7, three are relatively faint sources adjacent to 
bright sources and one is in a particularly crowded re- 
gion. The remaining three could potentially be pre-main- 
sequence objects. However, due to the lack of observed 
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Fig. 12. — Distance- and reddening-adjusted J — H vs. H — Ktor 
the G305.3+0.2 cluster region, statistically corrected for field star 
contamination. Error bars are comparable to the size of the points 
or smaller. Very few sources fall outside the region of reddened 
main-sequence stars by more than 2cr. 
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Fig. 13.— K vs H-K for the G353.4-0.36 cluster region. Note the 
clear separation in color between cluster and foreground sources. 
Sources Bl and B2 are denoted by filled symbols; Source B3 was 
not detected in H. 



gaseous emission from the cluster, we consider it unlikely 
that these are truly pre-main-sequence stars, and exclude 
them from the analysis along with the objects in the un- 
physical blue region of the color-color diagram as likely 
suffering from blending or a mismatch between sources in 
the different bandpasses. There are few enough sources 
in this region that we do not expect their inclusion or 
exclusion to greatly affect the IMF determination. 

3.2.2. G353. 4-0.36 

The J, H, and K color composite of G353.4-0.36 
is presented in Fig. [21 The youth of this cluster is 
immediately apparent from its heavily embedded nature 
and the dense molecular cloud that surrounds it. This 
region has long been known to be a site of massive 
star formation, and it has been studied extensively 
in the radio and sub-mm, includ ing continuum ob - 
servations at 1.5 GHz, 5 GHz i|Becker et alJ [1994), 
and 850 fim ijCarev et al.l 120001) as well as molec- 
ular line obse rvations in CS llGardner fc Whiteoakl 
EaSI, GO JWhiteoak. Otnmcek fc Renmel 119821 
HoGO llGardner fc Whiteoakl 11978ft . HNCO 
ijZinchenko. Henkel. fc Maa l2000|) (identified as a 
dense molecular core ), and SiO {Hariu et al. 1998). 
These signatures of ongoing star formation, combined 
with the strong nebular emission still present around 
the sources observed spectroscopically, suggest that 
the cluster is quite young, without main-sequence 
stars. Many of the continuum and molecular line 
observations quote slightly different positions for the 
source peak, and sources separated by several tens of 
arcseconds are all identified with the IRAS point source 
17271-3439. Since the beam sizes in many instances 
are comparable to the size of the NIR-bright nebulosity 
and to the separation between sources, it is likely that 
the extended source measurements are observing the 
same complex, which may peak at different locations 
in different w avelengths. Many of the radio da ta are 
tabulated by iGhan. Henning. fc Schreverl (^96), who 
identify a massive YSO in the region based on the 
IRAS colors. It is obvious from the NIR imaging that 



this source is not a single point source; in addition to 
the NIR source s, there are at least four separate sets 
of masers (e . g. iCaswell et all l2000l lArgon et all l2000l 
iVal'tts et al.l l2000[). one of which is associated with 
an UCHH (|Forster fc Caswelll UoOOi . Positions of the 
masers are indicated in Fig. [7| We note that the masers 
occur in regions which arc heavily extincted in the 
near-IR. OH, H2O, and CH3OH masers are all known 
in the region; the latter in particular are indicative of 
ongoing massive star formation. Clearly the sources 
visible in the near-infrared are only the tip of the 
iceberg, with other massive stars still in the process of 
formation. Higher-resolution maps at radio and sub-mm 
wavelengths are necessary to obtain a full understanding 
of this region. 

In the region of the large dark molecular cloud, only 
foreground stars are visible. This implies Ay > 50 in 
order to completely obscure the stars even in K, assum- 
ing a K-band detection limit of 17 and a distribution of 
K magnitudes similar to the rest of the field. The less 
heavily extincted region in which the cluster is visible 
in the near-IR must have been partially cleared out by 
stellar winds and ionization from massive stars. The rel- 
ative position of the NIR stars and the methanol masers 
(which lie in regions of higher extinction) suggest that 
we are observing stars nearer the main sequence which 
are emerging from the dust, while objects at an earlier 
evolutionary state are offset from this region, indicating 
ongoing star formation. 

The color-magnitude diagram of the G353.4-0.36 clus- 
ter is presented in Figure ITS! 

The cluster sequence is much narrower and more well- 
separated than in the G305.3+0.2 cluster, allowing for 
reliable separation of foreground objects based solely on 
H — K. Thus, we did not carry out a statistical removal 
of foreground objects for this cluster, instead considering 
only the objects well-separated from the foreground se- 
quence. Due to the high extinction toward this cluster, a 
large number of objects in the cluster area were detected 
only in K (shown as limits in Figure H~3*|) . The KLF is 
thus likely to be more reliable than the color-magnitude 
and color-color diagrams in determining cluster charac- 
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4. THE K BAND LUMINOSITY FUNCTION AND THE 
INITIAL MASS FUNCTION 

Onc e field stars have been rejected as described in 
< j3.2. II we can compute the KLF for both clusters. For 
the G305.3+0.2 cluster, which has more than 100 stars 
remaining, we additionally compute the initial mass func- 
tion (IMF) using two different techniques, the first using 
the KLF and the second using the color-magnitude di- 
agram and the spectroscopy of the massive stars. The 
KLF is commonly used to determine the IMF even when 
multi-color photometry is available; we take this oppor- 
tunity to test the robustness of this method and com- 
pare the results between this simple and commonly used 
method and the more involved method using the color- 
magnitude diagram. This will help to understand the 
uncertainties and systematic errors that may be a factor 
when only the KLF method can be used to derive an 
IMF. There were too few stars to robustly compute the 
IMF for the G353-0.4 cluster, so we compute only the 
KLF in this case. 

4.1. The G305. 3+0.2 Cluster 

To provide a robust determination of the KLF and the 
IMF, we must determine the completeness of our data, 
which we established by performing artificial star tests. 
Five artificial stars at a time were inserted into the clus- 
ter region; the small number was chosen to avoid sig- 
nificantly changing the crowding characteristics. IRAF- 
DAOPHOT was then run on the images to determine 
the number of artificial stars that were successfully re- 
covered. The procedure was repeated 50 times for each 
magnitude bin (Am = 0.5), for a total of 250 artificial 
stars added in each bin in H and in K . Figure IT~TI shows 
the results; completeness falls sharply to about 25% at 
H ~ 16.5, K ~ 15.5. We can compare these magni- 
tudes with the turnover in the "field luminosity func- 
tion" , which also probes incompleteness. The counts in 
the field turned over sharply at K ~ 16, in reasonable 
agreement with the artificial star estimate of incomplete- 
ness. 

4.1.1. The K Luminosity Function 

Knowing our incompleteness, we can calculate the KLF 
for the cluster. Figure El shows the uncorrected data, 
with the field "luminosity function" normalized to the 
same total number of stars overplotted for comparison. 
Figure 1161 shows the results after correcting for incom- 
pleteness by dividing the number of stars in each magni- 
tude bin by the recovered fraction of artificial stars. 

As expected, there is an overabundance of bright stars 
(K < 14.5) in the cluster region relative to the field. 
This is not an artifact of incompleteness; the complete- 
ness fraction at this magnitude is ~ 90%, and we expect 
incompleteness to be higher in the cluster than the field 
due to the effects of crowding. Using the number counts 
corrected for field star contamination (as discussed in 
M3.2.1I) and incompleteness, we fit a slope to the number 
counts in bins of AK = 0.5. We excluded sources fainter 
than K — 15.5 from the fit since errors in the incom- 
pleteness determination arc likely to dominate the num- 
ber counts. We derived a slope of 0.21 ± 0.06 for log N*. 
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Fig. 14. — Completeness fraction determined by artificial-star 
tests for the G305.3+0.2 cluster region (K = solid line, H = dashed 
line). 




10 11 12 13 14 15 16 17 



K 

Fig. 15. — K-band luminosity functions for G305.3+0.2 cluster 
(solid) and field (dashed), normalized to the same total number 
of stars. Note the peak is shifted to brighter magnitudes for the 
cluster. 
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Fig. 16. — Completeness-corrected K-band luminosity function 
for G305.3+0.2 cluster (dashed). The uncorrected KLF is over- 
plotted (solid) for comparison. 



This slope is somewhat flatter than the KLFs derived 
for more massiv e embedded clust e rs (e.g . 0.41 ± 0.02 for 
NGC 3576 from iFigueredo et all l|2002(l . 0.40 ± 0.03 for 
W42 from lBlum et all l)2000|) '). This suggests that this 
cluster is more weighted toward massive stars than the 
norm. 

4.1.2. The Initial Mass Function 

In order to better compare our results with the lit- 
erature, and to explore how much of a difference the 
use of multi-color photometry and spectra of the mas- 
sive stars make in the determination of the IMF, we 
used two methods to derive an IMF for the G305+00.2 
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cluster. For both methods we use a distance to the 
cluster of 4.0 kpc (as discussed in M3.1.1|) . The first 
IMF-determination method, which uses only the KLF, 
is commonly employed even whe n multi-color photome- 
try and spectra a re available (e.g. iFigueredo et alJl2002t 
iBlum et al.ll2000() . This method is simply a transforma- 
tion from K magnitude bins to mass bins. To make this 
transformation, we first correct the observed K for dis- 
tance and extinction as discu ssed in i J3.1.1l Usi ng the 
stellar evolutionary models of iMevnet fc Maederl l)2003f) 
for solar metallicity, we relate the mass for each track 
to an absolute K magnitude for a star on the ZAMS. 
We transf ormed Lhni to K using the bolom etric correc- 
tions from lVacca. Garm anv. fe Shul| II1996D for the early 
spectral types andfMalagnin i et alJ l)1986j) for later spec- 
tral types. We then us e the intrinsic V — K col ors from 
IBessell fe Brett! l|1988j) for A-M stars and from iWeenerl 
l)1994[) for O and B stars. Finally we interpolate linearly 
between the masses available on the evolutionary tracks 
to find the masses corresponding to our magnitude bins, 
and fit a power law to the resulting mass function. Our 
resulting IMF slope is Y = — 1.5 ± 0.3, excluding the two 
lowest-mass bins where incompleteness is significant. 

Our second method of determining the IMF made use 
of our multi-color photometry and spectra to estimate 
individual extinctions and masses for cluster members. 
Spectral typing of the brightest cluster stars allows their 
mass to be determined fairly accurately for a given stel- 
lar evolutionary model. For the models described above, 
the mass of an 06V star is approximately 40 M Q , that 
of a BOV star is 15 M©, and that of a B2V star is 8 
M©. Although spectra are not available for most of the 
cluster stars, their masses, as well as extinctions to the 
individual stars, can be estimated from the accurate rela- 
tive photometry. The presence of an O supergiant in the 
cluster suggests that, while the most massive stars have 
begun to evolve away from the main sequence, none have 
yet gone supernova, and less massive stars should still be 
on the zero-age main sequence. Therefore, with the ex- 
ception of the few most massive stars (for which we can 
estimate masses from their spectral types) the cluster 
stars should be scattered around the zero-age main se- 
quence (ZAMS) primarily by differential extinction and 
rather than the effects of stellar evolution. We can then 
use the same models and conversions from theoretical 
to observed quantities described for the KLF method, 
with additional transformations from T e j f to H — K 
usi ng intrinsic colo rs from from IBessell fc Brettl i|1988f) 
andlWegnerl dl994f> an d from T e ff to spe ctral type from 
iReoolust et alJ (|2004ft or lJohnsonl ©66). 

This transformation from theoretical to observed quan- 
tities allows us to place the ZAMS on our CMD. If the 
cluster is sufficiently young that we can neglect the ef- 
fects of stellar evolution, as discussed in the previous 
paragraph, we expect the ZAMS will lie in the middle 
of the distribution of cluster stars. The ZAMS derived 
from the evolutionary tracks o flMevnet fc Maederl l|2003f) 
is overplotted on the distance and extinction-corrected 
CMD in Figure ^2 A significant number of stars are 
bluer than the ZAMS on this plot. We interpret these as 
stars which are less extincted than those used to deter- 
mine the average cluster extinction and thus have been 
over-corrected by using the mean extinction. The scatter 
of stars around the ZAMS suggests that the extinction 



varies across the cluster region. To correct for this, we 
move the stars along the direction of the reddening vec- 
tor until they lie on the ZAMS. If the resulting extinction 
differs from the mean cluster value by more than Ay — 5 
for a given star, we exclude the star from the analysis, 
as it probably suffers from poor photometry. Examina- 
tion of the color image of the cluster region (Figure ^1 
suggests that the variation in internal extinction in this 
region is relatively small; no dust lanes or color varia- 
tions across the cluster are visible to the eye. The exact 
value selected for the cutoff is somewhat arbitrary, but 
does not greatly affect the results; most of the sources 
thus excluded have derived extinctions that differ from 
the median value by Ay = 10 or more. 

Using the positions of the extinction-corrected pho- 
tometry along the ZAMS, we are able to more accu- 
rately place stars in mass bins. The endpoints of the 
bins were determined by the masses for which theoreti- 
cal tracks are present in the models we used. In order 
to have an adequate number of stars in each bin we con- 
structed bins using alternate tracks for the endpoints, 
rather than every track. The analysis was repeated 
for three different metallicities (Z = 0.1 0.02 001) us- 
ing the evolutionary tracks o f iMowlavi et al.1 Jl998l Z 
= 0.1). ISchaller et alJ Jl99l: IMevnet fc Maederl (|2003l 
Z=0.02) and ISchaller et all l|1992l Z=0.001L For the 
solar-me tallicity case the high-ma ss points (M > 9M©) 
are from IMevnet fc M aederl i2003l) wh ile the lower-mass 
points are from Schal ler et al .1 l|l 992ft . The difference 
in K for the two solar-metallicity tracks is always less 
than 0.1 magnitudes for the masses where the two sets 
of tracks overlap and for most masses is less than 0.03 
magnitudes. The high metallicity model should be con- 
sidered only as a limiting case since such a high metal- 
licity is not expected. The use of such a wide range of 
metallicities allows us to estimate the importance of this 
parameter on the final IMF determination. 

Given these sets of mass bins, for each metallicity we 
determine the number of stars per unit logarithmic mass 
interval after correcting for completeness. We then fit 
a power law to the data. The two lowest-mass bins 
(M < 2Mq), where incompleteness was significant, were 
excluded from the fit; uncertainty in the completeness 
correction applied could significantly influence the results 
in these mass bins. The resulting completeness-corrected 
IMF for the cluster is plotted in Figure IT7I The solar- 
metallicity models yield an IMF slope T = -0.98 ± 0.2, 
where the quoted errors are only the formal fit errors 
and should be considered an underestimate. The low- 
metallicity tracks yield T = —1.01 ± 0.2 for the same 
distance, suggesting that the cluster IMF determination 
is insensitive to metallicity for solar and sub-solar values. 
The Z = 0.1 tracks give T = -0.88 ± 0.15. 

4.1.3. Comparison of the IMF Methods 

The IMF slopes we derive using these two methods 
are marginally consistent within the error bars: T — 
-1.5 ± 0.3 for the KLF method, and T = -0.98 ± 0.2 
for the CMD + spectroscopy method assuming solar 
metallicity. Comparing these results individually to the 
Salpeter slope would lead to different conclusions, how- 
ever. The KLF method produces a slope that is very 
close to the Salpeter value, while the slope from the 
CMD + spectroscopy method differs from Salpeter by 
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Log(Moss) 

Fig. 17.— The completeness-corrected IMF for the G305. 3+0.2 
cluster. The fitted values are the masses derived from the 
ISchaller et alj 119921) stellar evolutionary tracks with a distance 
to the cluster of 3.4 kpc. The plotted error bars are given by as- 
suming the error in the number of stars in a mass bin is equal to 
the square root of the number of stars in the bin. The fitted line 
has a slope of -0.96. 



about 2a. While this difference in slopes could arise 
purely from statistical uncertainty, various systematic ef- 
fects should cause the KLF-derived slope to be steeper 
than the CMD-derived slope, as we observe. If the more 
massive stars are preferentially located toward the cen- 
ter of the cluster, as expected due to mass segregation, 
and if the extinction is higher in the center of the clus- 
ter, the mean extinction used in the KLF determination 
would be systematically low for the more massive stars. 
This method would then underestimate the masses the 
highest mass stars, thus steepening the slope of the IMF. 
Evidence that this effect may be at work is provided by 
the six brightest cluster members, all of which lie red- 
ward of the Z AMS in Figure ^] while the fainter mem- 
bers are scattered more evenly. A difference in Ak (and 
thus Mk) of 0.2 corresponds to 1-2 subtypes for massive 
stars and thus to a difference in the derived mass of at 
least 2 M Q . 

However, mass segregation can only provide a partial 
explanation for the difference in the IMF slopes; the stars 
for which we obtained spectra are not in the very center 
of the cluster (since crowding in the 2MASS image used 
to select spectroscopic targets prevented us from select- 
ing targets in the cluster core). An additional possible 
source of systematic error in the KLF method relative to 
the CMD method lies in field star rejection. In addition 
to the statistical field star rejection described in E I3.2.1I 
which was done before any further analysis and thus ap- 
plies to both methods, the CMD method has color-based 
field star rejection. The CMD method can reject fore- 
ground objects, which due to lower extinction are bluer 
than cluster objects, as well as background objects which 
are redder than the cluster. The KLF method includes 
these objects, which tend to be fainter on average than 
the cluster stars (since they are either at a greater dis- 
tance or are low-mass foreground stars) and thus finds an 
artificially high number of low-mass stars. We find the 
use of K photometry alone to derive the IMF is likely to 
produce an overly steep IMF in regions with significant 




Fig. 18. — The raw (solid line) and completeness-corrected 
(dashed line) KLF for the G353.4-0.36 cluster. Only sources with 
H — K > 1, which fall in the cluster color sequence, have been 
included. 



field contamination or variable extinction. 

4.1.4. Comparison With Other Young Stellar Clusters 

Most studies of young star clusters have found an ini- 
tial mass function consistent with the Salpeter slope 
of r = —1.35, generally wi t h uncertainties of 0.1- 



0.2 (e..e. IPVuereHo et a,l .1 l200l IMassev fc Hu7rto)IT99 
Hillenbra nd &: Carpenter) 120001: lOkumura et alJ 
including the extremely m assive R136 cluster in the LMC 
l)Massev Hunt erlll998D . A review of the results is pro- 
vided in lMassevI <|2003). In the case of NGC 6611, re- 
analysis of the same data by different authors has pro- 
duced drama tically different results; a n IMF of — 1 . 1 ± . 1 
was found by Hille nbrand et all l)1993[) , while a reanalysis 
with different treatm ent of extinction produced — 0.7±0.2 
l)Massev et al.lll995l . suggesting that the systematic ef- 
fects are at work in IMF determinations that are at least 
as important as the statistical errors, as we see in this 
work. Slopes significantly flatter than Salpeter have been 
rep orted for the Arch es cluster near the Galactic Cen- 
ter l|Figer et al.lll999|) . though later work s uggests that 
this re sult is an artifact of mass segregation: iStolte et alJ 
( 200l found a very flat IMF in the core of the Arches 
Cluster with a steeper IMF at larger radii, with an overall 
slope consistent with a Salpeter value. The flatness we 
observe in both the KLF and the IMF for the G305+00.2 
cluster using the CMD + spectroscopy method may sim- 
ilarly be due to mass segregation. In addition to the ex- 
tinction effects mentioned previously, fainter stars in the 
outskirts of the cluster could be indistinguishable from 
the field star density (especially given the high field star 
density due to the location of the cluster in the Galac- 
tic plane) and not fall within the cluster boundaries we 
employ. 

4.2. The KLF for the G353-4-0.36 cluster 

Completeness tests were performed for the G353.4- 
0.36 cluster using artificial stars as discussed above, and 
the completeness-corrected KLF is plotted in Figure ED 
Since the cluster is significantly less crowded and faint 
cluster stars less common, our detections in this cluster 
are nearly complete in K, even though our detection limit 
is brighter than in the G305.3+0.2 cluster. The turnover 
at K = 15.5 appears to be genuine rather than an ar- 
tifact of completeness. Perhaps lower-mass stars in this 
cluster are still more deeply embedded in the gas and 
dust, and thus we observe only the massive objects. 

Due to the small number of stars detected in this clus- 
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ter (N = 25, only 7 of which were detected in H) and 
to the early evolutionary stage of the objects, we did not 
attempt to determine an IMF for this cluster or to place 
objects on the ZAMS. While the individual objects we 
observed in the G353.4-0.36 cluster were intriguing and 
worthy of further study, we cannot analyze the cluster as 
a whole because there are so few objects. 

This cluster is a very promising target for study at 
other wavelengths more suited than the NIR to the study 
of YSOs and even earlier stages of star formation; the 
methanol masers and likely presence of massive YSOs 
suggest that several stages of massive star formation can 
be studied in this region. 

5. SUMMARY 

We present NIR images and spectroscopy of two 
young stellar clusters near radio sources G353.4-0.36 and 
G305+00.2. Our if-band spectrum of the brightest clus- 
ter star in the G305+00.2 cluster show it to be an 051b- 
061b star. Although the range of luminosities of super- 
giants prevents us from determining an exact distance, 
this identification suggests a larger distanc e than radio 
dista nce to the nearby methanol masers l|Walsh et al.l 
11997ft of 3.3 kpc. We also obtained spectra of early two 
B stars in the cluster. There was no nebular emission 
present in the G305+00.2 cluster, though a ridge of neb- 
ular emission, coinciding with 8 /im emission and masers, 
is present ~ 1' away and may indicate sequential star 
formation, with the masers and gas indicating ongoing 
star formation and the cluster the result of earlier star 
formation. We computed the KLF and IMF of this clus- 
ter, and found them to be steeper than that reported for 
most young clusters (r = —0.98 ± 0.2 for the more reli- 
able CMD-based method) but generally consistent with 
the Salpeter value. We find that computing the IMF 
based only on a single color of photometry is prone to 
systematic errors when differential extinction and field- 
star contamination are significant. 

Our if-band spectra of two of the three stars we ob- 
served in the G353.4-0.36 cluster were featureless, while 
the other showed CO absorption, which is consistent ei- 
ther with a cool foreground giant or a YSO. The ab- 
solute magnitudes derived based on the distance to the 
radio sources are too bright for these objects to be solar- 
mass YSOs. None of the objects were detected in our J- 
band photometry, making identification as YSOs based 
on NIR excess impossible. They remain candidate mas- 
sive YSOs, and observations at other wavelengths are 
needed to make a positive identification. The images of 
this cluster showed a region with intense nebular emis- 
sion embedded in a very dark cloud where earlier stages 
of star formation are progressing. 
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TABLE 1 

Magnitudes and spectral line identifications of sources in G305.3+0.2 for which we 

OBTAINED SPECTRA. NUMBERS IN PARENTHESES IN THE PHOTOMETRY ARE THE ERRORS IN THE 

LEAST SIGNIFICANT DIGIT. 



Star 


RA 

13 h n m 


DEC 
-62° 


J 


Photometry 
H 


K 


Spectral Properties 
Species A (/im) EW (A) 


Spectral 
Type 


Al 


41 s .04 


32' 56"8 


11.75(1) 


10.39(3) a 


9.58(3) a 


Br7 
NIII 
CIV 


2.166 
2.116 
2.078 


-5.7±0.6 
-2.7±0.7 
>-0.8 


05V-06V 


A2 


33 s . 88 


33' 27".l 


12.310(1) 


11.02(3) a 


10.34(2) a 


Br7 
He I 


2.166 
2.112 


6.2±1.2 
0.7±0.2 


B0V-B1V 


A3 


39 s . 50 


33' 28".2 


14.063(4) 


12.646(4) 


11.97(2) 


Br7 


2.166 


5.9±1.3 


B2V-B3V 



a 2MASS magnitude 



TABLE 2 

Sample photometry for all stars in G305.3+0.2 cluster region. 
Magnitudes < 11.5 are taken from the 2MASS Point Source Catalog, 
since the IRIS2 images are saturated. Sources Al, A2, and A3 are listed 
first followed by the remaining sources. Field stars that were removed 
before deriving the luminosity and mass functions are included. 



RA (J2000) 


DEC (J2000) 


J 


AJ 


H 


AH 


K 


AK 


13:11:41.040 
13:11:33.877 
13:11:39.503 


-62:32:56.77 
-62:33:27.12 
-62:33:28.17 


11.751 
12.310 
14.063 


0.010 
0.001 
0.004 


10.394 
11.020 
12.646 


0.027 
0.027 
0.004 


9.575 

10.342 

11.969 


0.029 
0.023 
0.018 



13: 


11 


:37 


.680 


-62 


:.'«: 


09 


.60 


11. 


.949 


0, 


041 


10 


,776 





,047 


10 


.185 





,037 


13: 


11 


:36 


.286 


-62 


:33: 


13 


,30 


12 


,488 





,010 


11 


,519 





,001 


10 


.655 





037 


13: 


11 


Al 


.620 


-62 


:33: 


17 


,40 


12 


936 





,029 


11 


,650 





,038 


10 


.953 





034 


13: 


11 


Al 


.111 


-62 


:33: 


18 


,39 


16, 


309 


0, 


036 


14 


,514 





,018 


11 


.142 





002 


13: 


11 


:39 


.268 


-62 


:33: 


2d 


,85 


13 


,467 


0, 


003 


12 


,095 





,003 


11 


.494 





,018 


13: 


11 


:39 


.439 


-62 


:33: 


03 


,63 


13 


,218 


0, 


,003 


12 


,018 





,002 


11 


.524 





,018 


13: 


11 


:43 


.767 


-62 


:33: 


26 


.39 


16, 


,065 


0, 


104 


13 


,312 





,003 


11 


.594 





018 


13: 


11 


:40 


.045 


-62 


:33: 


18 


,89 


13, 


,549 


0, 


,006 


12 


,236 





,006 


11 


.596 





018 


13: 


11 


:38 


.141 


-62 


:33: 


13 


.66 


13 


,428 





,003 


12 


,270 





,003 


11 


.745 





018 


13: 


11: 


:39 


.493 


-62 


:33: 


10 


.25 


16, 


021 


0, 


129 


14, 


,851 





,116 


11 


.860 





004 


13: 


11 


:40 


.458 


-62 


:33: 


03 


.65 


15 


,989 


0, 


019 


14, 


,766 





,029 


11 


.920 





003 


13: 


11 


:40 


.021 


-62 


:33: 


07 


.26 


13. 


,884 


0, 


,004 


12 


,591 





,003 


11 


.970 





018 


13: 


11 


:40 


.992 


-62 


:33: 


07 


.86 


14 


,258 


0, 


,005 


12 


.845 





,005 


12 


.123 





,018 


13: 


11 


;36 


.748 


-62 


:33: 


11 


.14 


14. 


.115 


0, 


,010 


12 


,898 


0, 


,005 


12 


.153 





053 


13: 


11 


:34 


.747 


-62 


:33: 


2d 


.02 


14 


,067 





,010 


12 


,742 





,002 


12 


.311 





029 


13: 


11 


:3d 


.525 


-62 


:33: 


11 


.13 


14 


.362 





010 


12 


,969 





,004 


12 


.334 





044 


13: 


11 


:40 


.031 


-62 


:33: 


11 


.38 


14. 


315 


0, 


010 


13 


,098 





,014 


12 


.518 





018 


13: 


11 


:40 


.433 


-62 


:33: 


23 


.29 


17, 


,628 


0, 


,108 


15 


,615 





,049 


12 


.698 





005 


13: 


11 


:39 


.217 


-62 


:33: 


08 


.44 


14 


.799 


0, 


,010 


13 


,578 





,007 


12 


.789 





,164 


13: 


11 


:38 


.080 


-62 


:32: 


59 


,21 


16, 


,267 





017 


13 


,972 





,005 


12 


.872 





018 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 



TABLE 3 

Photometric data for the spectroscopic targets 
in the G353.4-0.36 cluster. 



ID 


RA (2000) 


DEC (2000) 


J 


H 


K 


Bl 


17:30:27.8 


-34:41:28.1 




14.14 


12.85 


B2 


17:30:27.9 


-34:41:34.7 




14.59 


13.47 


B3 


17:30:27.8 


-34:41:40 






14.38 
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